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\ Abstract. Until now the problem of counting Latin rectangles m x n has been 

• solved with an explicit formula for m = 2, 3 and 4 only. In the present paper an 

I explicit formula is provided for the calculation of the number of Latin rectangles for 

any order m. The results attained up to now become particular cases of this new 
formula. Furthermore, putting m = n, the number of Latin squares of order n can 
also be obtained in an explicit form. 
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0. Introduction 

^ I A Latin rectangle m x n is a matrix with n rows and n columns the elements 

of which are chosen in [n] = {1, . . . , n} so that two elements are never the same, 
neither on the same row nor on the same column. It is said that such a Latin 
. rectangle has order m. Prom the definition it follows that m < n and that each row 

5^ i of a Latin rectangle is a permutation of [n]. 

Furthermore it is clear that it is always possible to standardize the first row 
making it the same as the permutation 12 ... n. In such a case we say that we are 
dealing with a "reduced Latin rectangle" . 

If we call the number of Latin rectangles m x n with L(m, n) and the number 
of reduced Latin rectangles with the same dimension with K{m, n), it is clear that 
L(m, n) = nlK{m, n). 

The problem of counting Latin rectangles has engaged several generations of 
mathematicians, but the results reached up to now, as we will see later, are limited 
to certain special cases. 

With this paper we intend to finally supply an explicit formula for the calcu- 
lation of K{m, n) for any value of the order m. 

The partial results attained up to now will result special cases of such a formula. 
The result obviously also allows the calculation of the number of Latin squares 
of order n, which represent the special case of m x n Latin rectangles in which m 
takes on its maximum admissible value n. 
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1. A BRIEF SURVEY OF RESULTS 



A Latin rectangle consists of m permutations of [n] which, taken two by two, 
don't have fixed points. It is from this point of view that the problem was initially 
studied by Montmort, Euler and Lucas. 

It seems that the solution in the simplest case m = 2, known as "derangement 
problem" , can be found going back to Montmort 1713 [6]. It consists of the number 
Dn of the permutations of [n] without fixed points given by: 



and equivalent to K{2,n). 

In 1891 Lucas expounded the famous "menage problem" which consists of 
counting the ways of arranging n couples at a round table so that men and women 
alternate and no husband and wife are adjacent to one another. The problem, 
examined since 1878 by Tate, was also studied by Cayley and Muir, however no 
satisfactory results were reached. 

The solution to menage problem is equivalent to the enumeration of all the 
permutations of [n] which are discordant with both the permutations 12 ... n and 
23. . .nl. The generalisation of the above mentioned problem — known as "the 
cyclical Touchard problem of index m" or "problem of m-discordant permutations" 
— proposes the counting of all the permutations a of [n] so that: \/ i & [n], a{i) — i ^ 
0, 1, . . . , m — 1 (mod n). 

It is also said that this problem is equivalent to the enumeration of the "very 
reduced Latin rectangles" of order m + 1, the number of which is indicated with 
V{m,n), that have the first m permutations aj, with j G [m], in the canonical 
form: aj{i) — i + j — 1 (mod n). 

The solution in the simplest case m = 2 was found by Touchard in 1934 [15] 
and consists of Touchard's famous numbers Un, equivalent to F(2,n), expressed 



For the subsequent case V{3, n) recursive algorithms have been obtained by Riordan 
[12] and by Yamamoto [18]. However an explicit formula was only provided in 1967 
by Moser [7]. In the case y(4,n), the biggest yet dealt with, there is only one 
recursive result by Whitehead in 1979 [16] and possibly an explicit formula by 
Nechvatal [8] also in 1979. 

However let us return to the more general and more complex problem of the 
calculation of K{m, n), which represents the aim of this paper. 

The first attempts at the calculation of K{3,n) go back to Jacob and to Ker- 
awala who, in 1941 [5], found a recursive formula. The following tidy explicit 



(1.1) 




by: 



(1.2) 
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formula for K{3, n) is, on the other hand, attributed to Yamamoto (see [1] and 



Furthermore, in 1944, Riordan obtained an expression of i^'(3, n) in terms of Touchard's 
numbers and subsequently, in 1946 [11], the well known formula: 



(with Uq = 1) which expresses K{?>, n) in terms of Dk and Uk- 

It is necessary to say that until now, in this line of research, no other progress 
has been achieved since, for m > 3, it has not been possible to obtain K{m, n) in 
terms of K{i,n) and V{j,n) with i,j < m. 

The case of i^(4, n), which is the most complex yet to be dealt with successfully, 
was only solved with an explicit formula in 1979, this was achieved independently 
by Nechvatal [8] and by Athreya, Pranesachar and Singhi [1]. 

Subsequently, in 1980, Pranesachar [10] and Nechvatal [9], by different means, 
found a way to express K{m, n) for any value of m by means of the Mobius function 
of the lattice of partitions of a set. The limit of these research works is that they 
take the calculation of K{m, n) back to the enumeration of other combinatorial 
objects, such as the partitions of an integer, for which no explicit formulas are 
known, and thus don't allow an explicit formula for K{m, n) to be obtained. A 
further tidy result of this type was achieved by Gessel in 1987 [3]. 

Until now, then, no explicit formula is known which permits the calculation of 
K^m, n) whatever the value of m. 

We would like to conclude this section remembering that another interesting 
line of research tried to get asymptotic expressions of K{m, n). The first significant 
paper of this kind is attributed to Erdos and Kaplansky [2] in 1946, subsequent 
results were obtained by Yamamoto in 1951 [17] and by Stein in 1978 [14]. 

Finally, in recent years, Godsil and McKay [4] achieved an asymptotic valuation 
of y(m, n). 



Very useful concepts in the study of permutations without fixed points are 
those of board and of rook polynomial. 

A board is a nonempty subset of P x P (P = set of positive integers), the 
elements of board are called squares. Considering a board C it is usually indicated 
with rk{C) the number of different ways of placing k non-attacking rooks on it. 
The rook polynomial of C in the symbolic variable x, that we'll call R{C), is given 



[10]): 



(1.3) 




(1.4) 
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If we write &n for the set of all permutations of [n], then every a e S„ can 
be thought of as a board, called the "graph" of a, the squares of which are the 
couples {i,a{i)) G [n]. It is furthermore obvious that m permutations, two by 
two without fixed points, make up a board of m ■ n squares, if in such a board 
the first permutation is the identical one a{i) = 1, it will hereon be indicated with 

We furthermore state, to have a greater number of symbols, that a number 
put up to the right of a symbol doesn't denote a raising to a power of the same, 
but it acts as a new symbol (therefore it isn't the square of a). When we want 
to indicate a raised to m, we write (a)"^. 

We will call the board, included in C^'^\ formed by the g-ih permutation 
of C("^) {g e H). Since, as has been mentioned previously, to speak about C'^"^-' 
it is the same as to speak about a reduced Latin rectangle of order m, we will refer 
often to as the ^'-th line of C^"^^ (which is not to be confused with the g-th row 
or column of C^"*^ like a board that are different things). We can even say that a 
subset of has "grade" g. 

At this point let us remember a classic result which joins the rook polynomials 
to the permutations without fixed points. We consider the permutations as boards 
and, taking the board 5 C [n] x [n], we indicate with Ng{B) the number of permu- 
tations of [n] which have exactly s squares in common with B. And so giving us 
the following tidy relation: 

s ^ ' 

for the proof of this see [13] chapter 2.3. 

Now let us introduce some conventions in the use of symbols. The sets are 
always indicated with capital letters and the number of the elements which make 
them up with the corresponding lower case letter (therefore a = \A\). C{A) will be 
the complementary of the set A in the universe set. 

j'ini) ^^Yl represent a generic system of independent rooks — that is which 
don't attack each other — put on C^"^) and T^, with g E [m], will be the part of 
the system contained on the g-th line of C^"*) (T^ = T^"^) n C») and thus it will 
be: T(™) = U T2 U ■ ■ ■ U T"^. 

Furthermore, if A C C^'^\ we will say that TZ^A) is the projection for rows of 
A in C^'^\ consisting of all the squares of C^"^) which have any square of A in their 
own row (obviously A C TZ{A)). We will also say that IZg^A) is the projection for 
rows of A on the ^-th line of C^"*) and we wiU put ng{A) = 7^(A) n C^, with the 
consequence that: 'R-{A) = 7li{A) U • • • U TZmiA). Similarly, speaking of columns 
instead of rows, we can define C{A) and Cg{A). 

Finally we define the set X(A) = TZiA) r\C{A) as the "impression" of A and the 
set 0{A) = n{A)UC{A) as the "shadow" of A. Thus 0(T('^)) wiU be the set of aU 
the squares of C^'^^ subject to the attack of any rook of T^"^^ and which therefore 
can't contain other independent rooks from those of T^"^) . 
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To indicate the number of ways in which the set A can generally be arranged, 
considering the restrictions which have been imposed on it, we will write tt{A). So 
we shall obtain that K{m,n) = 7r(C*^'^^). 

As it is known, a partition in k blocks of a set A is formed by a collection of 

k 

non empty sets Aj, with i e [k], two by two disjoint and such that (J^ Ai = A. We 

1 

will indicate with n(^) the set of the partitions of A; if tt e IK^) tt has k 
blocks, we say that |7r| = k; finally we put Yin — 11 (N)- 

Let us also remember that the refinement of two partitions tti and tt2, with 
TTi, 712 G IK^)' partition of A consisting of all the non empty intersections of 

some block of tti with some block of 7r2. 

If X = {xi, . . . , Xs} is a set of variables, we put, for economy of space: 

("■) = ( " ) 

and furthermore: 

^ X, ; n^= n ; n^'= n 

Xi£X Xi£X Xi€X 

and, ranging each Xi within its own domain: 

E = E-E- 

We will indicate with {n)k = n{n — 1) . . . (n — /c + 1) the falling factorial of n and 
with (n)fc = n{n + 1) . . . (n + /c — 1) the raising factorial of n. Now, given that 
{n)n — nl, we intend to put {n)n — n\ even if in other literature this symbolism has 
been used to indicate the subfactorial of n. 

Let us conclude this section with some recalls relative to the permutations of 

If o" e ©n and a{i) — ai \/i & [n], we can also say that a corresponds to the 
word ai . . . an- It is well known, see [13] pages 17 and following, that a can be 
shared in an unambiguous way in the product of disjoint cycles on the elements of 
[n] and that it can have a "standard representation", which we will indicate with 
Si . . . Sn, writing a) the elements of each cycle with the largest element first and 
b) arranging the cycles in increasing order of their largest element. In such a case 
each cycle will start with a left-to-right maximum i. e. with an element Si so that 
Si > Sj for each j < i. 

Another way to describe a can be achieved by indicating with bi the number 
of elements j of its standard representation on the left of i with j > i and defining 
I (a) = (bi, . . . , bn) the "inversion table" of a. In fact it can be easily proved, see 
[13] Proposition 1.3.9, that between the a e and the I {a) there is a bijection 
and furthermore that: < 6j < n — z, Vz e [n]. 
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3. The associated partitions 

In this section we want to show how the computation of K{m, n) can be taken 
back to the enumeration of a double system of partitions of [n] . 

The first step in the argument is to reiterate the method of calculation by- 
means of systems of independent rooks expressed by the formula (2.1). 

Let us suppose that we want to determine the number K{m, n) of all the 
possible C*^"^-* and that we have already counted all the possible arrangements of 
the first m — 1 lines C^'^~^\ the number of ways in which C"^ can be arranged 
can be obtained easily, by means of the formula (2.1), once rk (C*^"^"-*^)), for k = 
0, . . . , m — 1, are known. 

In fact, if r|^^| = Ty^_^ U • • • U T[^~^j^ is the independent generic system 

of rooks on C^"^~^^ and therefore t^^^]^i — ti/m-i ' — '"^i/m-i ~ ^' have 
that: 

(3.1) n{C-) = (-1) (- - A7r:^-^y■- (^iM-i) • 



Now we are trying to calculate 7r(C"^~^) with the same assumptions. Now 7r(C"^~^) = 
7r(T.";-i . )7r ( C"^-i - T:^-^ , ] and, if we consider the generic T^T'^l Cl( TV"-^ , \ 

and T^T'^l C I fc"^-i - TT]-^ .V we wiU have that: 

2/m—2 — \ 1/m—l / ' 

(3.2) n{C-')= Y: E (-1)*^^"^-H'v--i - ^V^^^^^ 

,(m-2) ,(m-2) 
l/m-2 2/m-2 

\ l/m—2 } ^ ' \ 1/m—l 2/m— 2 I \ 2/m—2 I 

Repeating the argument for the subsequent line C"^~^, we see that: 

(3.3) = TY".-^ „ U ( TV",-^ , - T„";-2 Tm-2 ^ ym-2 \ 
V ' l/m—2 \ 1/m—l 2/m— 2 I I 1/m—l 2/m— 2 I 

I I / rpm— 2 ry-tm- 2 \ I I ( /~tm—2 /rpm—2 , , rpm— 2 , , rpm — 2 \\ 
^ l^^2/m-2 ~ ^1/m-l ) ^ ~ \^l/m-2 ^ ^1/m-l ^ -'2/m-2^ ) 

and, calling: T^^'X T^T'X T^T'X T^T~% T^T'^'l the generic T^^-^) in- 

■ o 1/m — 3' 2/m — 6^ i/m — i^ A/m — 6^ o/m — 6 " 

eluded in the impression of each of the five components of the union of which at 
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(3.3), we will be able to say, referring always to (2.1), that: 
(3.4) 

_ / rrim — 2 ^ T-im— 2 I „ / rpm — 2 rpm — 2 \ „ / /^m— 2 rr\m — 2 

' 1/m— 1 2/m— 2 / \^ 2/m—2 1/m — l J \ 1/m—l 

5 

. - , * i/m — 3 

U T™-^ . U T^"^ 



■l/m-2 ^ 2/m-2 



)= """'(C/„i.-C-'3)' 



^(m-3) 
i/m — S 



j \rpm-2 _ rpm—2 i _ ,{m-3) \ | / |^m-2 p, rpm-2 i _ ,(m-3) \ | 
" y l/m-l 2/m-2l '^2/m-3y ' y 1/m-l ' ' 2/m-2 1 '^3/m-3j- 

( \rpm—2 rpm—2 i ,(m— 3) \ i n^m—2 / rpm — 2 , , rpm—2 i i rpm — 2 \ i 

■ \^l^2/m-2 ^1/m-ll ''4/m-3)-\\^ y- \ j m-\^ ^ \ j m-2^ ^ 2 j m-2) \ 

_ Arn-3) \ , / rp{m-3) \ f rp(m-3) \ f rp{m-3) \ ( rp(m-3) \ f rp{m-3) \ 
^hlm-3)-^ y-\lm-3 J ^ l^^2/m-3 J ^ \-^3/m-3 J ^ l^^4/m-3 J ^ \-^5/m-3 J ' 

Continuing in this way in order to count all the possible arrangements of the line I, 
we should take all the possible subsets of rooks which are situated on it, considering 
the partition refinement of consisting of all their non empty intersections Tj, 

pi 

j = 1, . . . ,pi, and of the complementary of their union Tq = — Uj^j finally 
choosing pi + 1 systems of independent rooks T^^-i' J — ^i^j ■ ■ ■ ^Ph with ^^7-1 — 

x(t;). 

Then we will have: 



Z^J j/l-l 



(3.5) .(C) = UMTi) = E (-1) ' n, (4 - *?/;-l)'-(2;';r-;) 



and we can conclude that: 
(3.6) 



^(i-i) 



m Pi 



j/l-l m pi 



K(m,n) = n. AO = E (-1) ' ° n. n, i!', - >^j^^'-'^fii^ 



M-i) 2 



The situation, therefore, seems to be somewhat complex, but an idea which al- 
lows us to control it is that of identifying any subset of C*^"^-* by means of its two 
projections, for rows and for columns, on the main diagonal = i G [n]}. 

Given a set of grade g , we consider in fact TZi{A^) and Ci(A^). Now, if 
g = 1-1 TZi{A^) = Ci{A^) = A^ and projections and set coincide. If, on the other 
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hand, ^ > 1, then the set determines TZi{A^) and Ci{A^) in one way only. 
Viceversa given Af , Af C with \Af\ = \A?, \ = \A3\ = p, A^ will be one of the p\ 
permutations of the square board TZ {A\) fl C (A|). Furthermore, if i? is a board of 
forbidden positions for A^ and we set i? = S Pi 7^ (Af ) n C we will have that: 

T:{A3) = j2k{-lf{v-k)\ru{B) . 



In the light of this new approach, a system of independent rooks T*^') = U 
••• U determines the subsets W = n{T') and = C(T*), for i G [/], and: 

= n{C{T^^^)) = c([j,R^ and (7° = C (C(r«)) = C ([jiC^y Thus it 

characterizes two partitions, each one with I + 1 blocks, the {R^} and the {C*}, 
with i = 0, . . . , Z, the first for the set of the rows and the second for that of the 
columns. 

Now, if we return to the computation performed with (3.5) of all the possible 
arrangements of C*^'-*, we will have that ^j/;^] characterizes a partition in / blocks 
Rj/i-i = ^(^j7z-i), with i = 0,...,/ - 1, of the block 7^(^j) of the partition 
{TZ{Tj)}, j = 0, . . . ,pi, in pi + 1 blocks of the set of the rows. 

Consequently, if we put together all the elements of grade i of the various 

partitions R^- n_i and we set R\ = IJji?* we will obtain that {i?;}, with i = 

0, ...,/ — 1, is a partition in / blocks of the set of the rows of our board. Similarly 
we can construct the partition {C|}, with i = 0, — 1, of the set of columns. 
Repeating this for each line from the m-th to the second, we will eventually have 
m — 1 couples of partitions of the set of the rows and of that of the columns: {-Rj} 
and {C*]}, with j = 2, . . . , m and i = 0, 1, . . . j — 1. 

Intersecting each of these partitions with we obtain as many partitions of 
= with i G [n]. In such partitions the projections for rows and for 

columns of the sets of grade 1, which belong to C^, obviously coincide. 

Calculating the number of all the possible arrangements of these 2(m — 1) 
partitions, respecting the condition that the projections of the sets of grade 1 
must coincide, is equivalent, from what has been said, to calculating the product 

rii rTi^r \ T\n_[ ) which appears in (3.6). 
2 ^ ^ 

To do this it is natural to consider the partition refinement of the m — 1 
{Rj n C^} partitions and that of the m — 1 {C] fl C^} partitions. Putting by 
analogy R^ = Ci — , the blocks of the refinement partitions will be given from: 
^a,„.....ai = R^"' n R'^Si n---nR^\ with O < < j - l \/j e [m], and from 
Cf3^,...^/3j^ — n C^Si n • • • n Cf^ with similar limitations on the indices (3j. 

We will then have: 

m pi j-1 j-1 

2 
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We will also say that the pair of partitions {Ram, -;cei} and {(^^^...../Si } of [n] 
is "associated" with the collection of systems of independent rooks T^'^-i' with 
I = 2, . . . ,m and j = 0, . . . ,pi. 

4. The blocks of the associated partitions 

Before being able to develop the calculation of (3.6) using (3.7), we must exam- 
ine closely the meaning of the indices am, . . . , cti and jSm, • • • , /5i which respectively 
mark the blocks Ram,...,ai and C^^^...^^^ . 

In order to do this, we first of all define the concept of "covering". Taking 
the index a^, with ctj > 0, we say that ai "covers" ckq^ and we write ai h ckq^ or 

K{ai) = tta,. 

/,From the definition it immediately follows that: 

a) if cti = it doesn't cover any other index, since ao doesn't exist; 

b) if Q!j h as then s < i. 

Thus, if CKs > , applying to it more times the function of covering re, you will 
always arrive at an index of value 0. 

On the contrary, taking an index ai of value 0, we can consider all the indices 
as which have the property ag H ai (that is K~^{ai)). Repeating this procedure 
more times, we get all the indices which, with a finite number of applications of the 
function re, finish in ai. 

If we suppose that ah = and we put K^{as) — as, we will be able to define 
Zh = {aj I = and 3 G N so that K^{aj) = a^} which we will call the "com- 
ponent h" of the indices am, ■ ■ ■ ,ai. Furthermore, as ai — 0, we will have that 

In this way, we obtain a partition of the set of indices {am, ■ . ■ ,ai\ in blocks 
made up from Z/^, with h e [m]. 

We will say that such a subdivision represents the structure of the indices 
am, ■ ■ ■ ,ai and we will write that cr(«m, . . . , ai) = ZiUZ^^^- ■ 'UZ^^ with 1 < Z2 < 
■ ■ ■ < Za < m. We will also write ai{am, ■ ■ ■ , oi) = Zi and C(ctm, ■ ■ ■ , cti) =number 
of the aj which are equal to zero. It is clear that, if Z = {ai^ , . . . , ai^} is a 
component, then ctj^ = 0. 

Moreover we will put, to be brief Ra^ = Ram,...,ai, C/3. = Cp^^,,,^^^, aj = 
am, ...,ai and Pj = (3m, 

It is necessary to pay attention to the fact that Zh is not only a subset / C [m] , 
but a subset of the indices aj, for j e /, each with its own value. 

The following result allows to count the number of the Raj at the base of the 
structure of their indices. 

4.1 Proposition. Let (j{am, ■ ■ ■ , cii) = ZiU Z and Z = Z^^ U ■ ■ ■ U Z^^, then: 

a) if we suppose Z to be variable, the number of the possible sets of indices will 
be: r;')z\. 

b) if, on the other hand, we keep Z constant, then the possible aj will be (m — 
l-z)\. 
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In fact, to determine Z we will have, first of all, have to choose the z places of its 
indices in the set {2, . . . ,m} and this can be done in C^"^) ways. Furthermore, if 
the selected indices are Ofj^ , . . . , (with ji > ■ ■ ■ > jz), it can be seen that aj^ 
must be equal to 0, aj^_-^ can assume the values and aj^ and so on. Therefore 
the last index has only one possible value, the penultimate two values etc., thus all 
the possible ways to attribute a value to aj^ , ■ ■ ■ , ocj^ are 1 - 2 ■ . . . ■ z = z\. And this 
proves a). 

If, on the other hand, Z is fixed, the places of the indices of Zi arc also fixed. 
Now the last index of Zi on the right ai can only have the value 0, the penultimate 
only 1 and thus, for an argument identical to the previous, the number of possible 
values of the m — z indices is equal to 1 • 1 • 2 • . . . • (m — 1 — 2;) = (m — 1 — z)\. And 
thus b) too is proved. 

It is also possible to calculate the number of possible aj , in terms of the data 
of singular components with the following result, which we shall just state. 

4.2 Proposition. If a{am, ■ ■ ■ , ai) = ZiU Zz^^- ■ Z^^, then\Zi\ and\Zzi\, with 
i = 2, .... a, constitute a partition of the integer m in which the number of parts 
equal to s will be Xg. The number of possible aj with this structure will therefore 
be the same as: 



1^12^2 ...m^- Ai!A2!...A^! 



So, if as \- cxij we have that is a subset of TZi ^|J^ '^x/i-ij^ with x that ranges 

in a subset of {0,1, ... ,pi}, and thus it lies in the projection for rows on the first 
line of a set of rooks of grade included in the impression of a set of rooks of 
grade ag = I. If instead ag doesn't cover ai, then the set of rooks of grade ck/ is 
not included in the impression of the set of rooks of grade ctg and so the number of 
elements in their intersection varies according to the variation of the projection for 
rows or for columns. 

From this, it follows that if, using the symbolism of section 3, we take TZi (Tj) — 
i-i 

T^ii U we see that it will be composed of the union of all the i?Q, . with the 

1 ' 

same Zi component. Viceversa, if we fix the Zi component and make the other ctj 
vary in all possible ways, we obtain a collection of sets R^- the union of which will 

be equal to 7^l (Tj) - 7^l(^jJ7iJ) for some j. Furthermore, if / = 1, since ^y^^] 

doesn't exist, the union of all the Ra.^ with the same Zi will be given by TZiiTj) 
for some j. 

Naturally the same argument is true for the sets Cp. and the components of 
the indices (3j. 
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5. The enumeration of Latin rectangles 

Now let us try, applying the contents of the previous section, to give an explicit 
form to (3.6) in terms of the data of the two associated partitions {Raj } and {Cj3.}, 
and that is in terms of the sets of variables R — {vaj} and C = {cp^}. 

First of all, we observe that, if i > 1, - tf~]_\ = |Tj| - |Tjj7_^J| = |7^l(Tj)| - 
\T^iiT-^/J-l)\ = l^i(ri) -^i(^j)7_^l)l but, following what was said before, 7^l (Tj) - 
7?.i(Ty7^|) is formed, in such a case, from the union of all the sets Raj with the 
same Zi component and viceversa. 

Therefore, if we put, V/ e [m], Q{Zi) = {vaj \ (^i{a.j) = Zi}, Q{Zi) = {cp. \ rp. e 
Q{Zi)} and q{Zi) = ^ Q{Zi), thus we have that, if i > 1, j G {0, . . . ,pi} exists so 
that: 

(5.1) ?(zo = 4-t5'/7!l 

and viceversa. 

Now, we will compute the product Yli Ylj "^(^ivz-i) that, for (3.7), is the same 

2 ■^^ 

as n\. n\.7r(i?a,)7r(C^J . 



The first partition {Raj } can be chosen in a completely arbitrary way and thus 
the number of its possible arrangements is given by the multinomial coefficient (^) . 
The second partition is, on the other hand, subject to some restrictions. 

First of all, for I > 1, T^^'^^ C T(T]) and so \C^{Tj) - Ci{T^'-}l)\ = |7ei (Tj) - 

TZi{Tj''^~y^)\ = q{Zi) and since, following the same reasoning as we have already 
done, iCi (Tj) - Ci(Tjj7_'J)| = E QiZi) we have that: 

(5.2) Yl ^(^') = ^(^') • 

Furthermore, taking a generic set of rooks of grade I it is clear that |Ci(T')| = 
\n^{T% 

Now, if / = 1, then C Ci and even Ci(T') = 7li{T)\ but, for the reasons 
stated in section 4, 7li{T)^ is made up of the union of all the R^ with the same 
Zi, and the same argument is valid for Ci(T^), therefore: 

(5.3) U Ra,= U C^, 

ai{aj)=Zi ai{(3j)=Zi 

and: 

(5.4) Yl ^(^i) = ^(^i) • 
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Furthermore, the restrictions (5.2) and (5.4) imposed on cp. imply that, for Tj, 
withi = 0,...,p«, |Cl(Tj)| = |7^l(Tj)|. 

This can be easily proved for complete induction on I considering that, if / = 1, 
the result has already been expressed by (5.4), while, if Z > 1 and we suppose that 
we have already proved this Vj e — 1], it follows from the consideration that: 

irjl = Ei I Tj/^-il + 1^1 (^i) - ^1(^77-1)1- 

Therefore there are no other restrictions on c^^. , apart from those expressed by 

(5.2) and (5.4). 

If we now group the Cp. sets on the basis of the value of their component Zi, 

(5.3) allows us to state that: 



on the condition, however, that the C variables also respect the restrictions imposed 
by (5.2). 

Let us finally examine the tj/i^^i which appear in (3.6) as exponents of —1. 

Pi 

Now, for (5.1), if i > 1, t^l-]_{ = t] - q{Zi) for any Zi and so E,*-/7-i = 

pi 

Yli3^\ ~ Yliz ^i^i) — ~ E^ q{^i)- Therefore, being I = 2, . . . ,m, the exponent of 

m 

— 1 in (3.6) will be the same as n(m — 1) — Ej E li^i)- Furthermore, since, as 
we have already seen, n = J2 ^ ~ ^(^i); adding and subtracting n it can be 



expressed by: nm - Y.i Qi^i)- 

1 ' 

Finally, set W = {r^j \ C(cij) odd} and considering that every Tq,^ variable 

m 

compares C(<^j) times in E( E q{^i)j we will have that the exponent of —1 can 

1 ' 

be substituted by nm — E ^) since the even multiples of Tq,^. can obviously be 
omitted. 

Using all these results in (3.6), we obtain the following remarkable result: 
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5.1 Theorem. 



nm-\-y i> q(Zi) 

(5.6) K(m,n) = E c M) ' (jn^II^^^^)' 



m 

III Viz, i-QiZi)) i 

-(_-[\n{m-l) Y" V -5 

=(-ir"E.E»(;)n.n.M)'<^"(g,'))-n 

1 C 



,C(/9,)-l 



Where, by analogy with the preceding symbohsm, we have set Zq — since 
ao doesn't exist. Thus q{ZQ) = ^ = n since the elements of Q{Zo), not being 
subject to any restrictions, are all the elements of R. 

So (5.6) is an explicit formula for the computation of K{m, n) in 2m! variables 
R and C, while q{Zi) with I e [m] and W are sums of particular subsets of R. 

This therefore represents the result which we proposed to achieve with the 
present paper. 

The C variables, in contrast to the R variables, are not, however, between 
their independent since they must be subject to the restrictions '^Q{Zi) — q{Zi) 
for I e [m]. 

If we want to limit ourselves to considering only independent variables we can 
proceed as follows. 

For each Zi component we indicate with d{Zi) the variable with (7i[(3j) 
= Zi and all the indices (3j which are different from those of Zi equal to zero, and 
we put D = {d{Zi)}. 

Now d{Zi) = q{Zi) — J2 (^Pj thus the variables of D can be obtained 



=/3j eQ(Zi)-D 



from those of C — D. 

Furthermore, if cri{Pj) = Zi, then cp^ e Q{Zi) and so C/j^. < q{Zi). Thus, if we 
put iif3. = min^^(^^^_^g {q{ai{pj))), then ^ cp. eC - D,we wiU have that 0/3. < up. 
and such a restriction guarantees that d{Zi) > 0. 

Using this new symbolism (5.6) can be rewritten like this: 

m 

(5.7) A'(„..n) = E,.,E,,(-i) ^ U) n(g-i>)!ng! 



14 



AURELIO DE GENNARO 



with Taj e R and cp. E C — D. 

Now if \Zi\ = s, for the Proposition 4.1, the possible d{Zi) are {^){s — 1)!. 
Furthermore, if s = 1, all the d{Zi) coincide with the cp^ which has all the indices 
at and so, in such a case, instead of ('7)0! = m we only have one distinct element 

m 

and |L»| = ^4-)(s-l)!-(m-l). 
1 

Thus, in (5.7), other than the ml independent variables R, there are the ml + 

m 

m — 1 — C^) (s — 1)! independent variables C — D. 
1 

6. Simplifications of the formula 

We have seen that (5.7) needs 2m! + m — 1 — Yls{^)is — 1)1 independent 

1 

variables for the computation of K{m, n). 

It is possible, though, to effect two types of elimination among these parameters 
which allow us to reduce their number considerably, even though this fact makes 
(5.7) lose its symmetry. This is obviously important when we would like to calculate 
concretely K{m, n) for m and n prefixed. 

Let us therefore examine the two possible reductions of the independent vari- 
ables R and C — D. 

A) We consider r^^ and Cq,^. with C(Q^i) = 1 and so with cr(aj) = Zi. In such an 
assumption q {a{aj)) contains a unique element and so, for (5.4), c^. = r^. and, in 
(5.7), Cq^. ! is simplified with q {a{aj))l = Vajl- As far as Tq,^. is concerned instead, if 
we put: Fo = {r^^. | Ci^j) = 1} , fo = J2Fo,Qo = R- Fq and Fq = {c«.|ra^. e Fq}, 
we will have that the variables of Fq don't appear in any set Q{Zi) with I > 1 and 
that: 

?(^)=?|:(.o:^J^?|:(;o)(i:)("«/°) 

" Qo Qo 

= yy((m-i)!)^°— ^ 

^V'' /o!nQo! 

Qo 

since, for the Proposition 4.1, |i^o| = (^ ~ !)'• Furthermore the Fq appear among 
the exponents of —1 with their total /q. The 2(m — 1)! variables of Fq and of Fq 
can therefore be substituted by /q. 

B) Let us now consider the Tq,^. and Cq,^. with (j{aj) = ZhU Zi (and so ({(^j) = 2) 
and mm{zi, Zh) = 1 and put, V.s' G [m]: Fg = {vctj \ C{ctj) = 2 and |a"s(aj)| = 1}, 

fs = E Fs, Q{Zs) = Q{Zs) - Fs, Fs = I r«^. G F,} and q(Z,) = g(Z,) - 
First we observe that, if Tq,^ G Fg, it doesn't appear among the exponents of — 1 
since Ci<^j) is even. Now, if cr{aj) = ZgU Zy, \Zy\ = m — 1 and so q{Zy) has only 
one element and, for (5.2) and (5.4), = ray Therefore in (5.7), if Cq^. G -Fg, Cq,^. ! 
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is simplified with q{Zy)\ = Ta^^-- Moreover, Vs e [m]: 

p Is- \^sj Js- 

since |Fs| = (m — 2)!, and so also the Fs and the Fs are eliminated and substituted 
by fs- We must, however, by careful because, if m = 2, Fi and F2 are equal. 

m 

In conclusion, putting R = R—[JiFi, C = {cq, | Tq, E R} and D — {d {ai (aj)) | Cq, . e 
C} we have that (5.7) transforms itself into: 
(6.3) 

K{m,n) = J2J2J2nJ2c-ni-^) ((m-l)!)^°((m-2)!)i 

fo fs 

It is, however, possible to accomplish a further step to simplify (6.3). In fact, 

m 

putting / = ^s/s, we have that: 
1 

(6.4) 

Jl Jm J- /I /m J- 

m 

// + Es^(^.) + ^-i^ 



and so (6.3) becomes: 



n n n I^Pj 



/oE/E7.E^-^(-i) ((^-1)!/° 



((m-2)!)^ 





n 



jo,f,Rj n(c-^)!n^! 

if + J2sQ{Zs)+m-ly. m m-2 m 

(EM+.n-i)/ n-n,.,^(^o!n.^(^.)! . 
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From the independent variables R we have so eliminated the (m — 1)! of and the 
(m — 2)! of each Fg, with s e [m], and therefore = ml — (m — 1)! — m{m — 2)! = 
m! - (2m- l)(m-2)!. 

The C variables have undergone the same reduction. However it is necessary 
to add /o and / and subtract the D, which are as many as the components Zi 

m-2 

with 1 < \Zi\ < m — 1, and so equal to ^ hi"h)i^ ~ 1)' P^^^ '^^j with all 

2 

the ttj indices equal to zero (which is determined by the m equivalent restrictions 
q{Zs) = E (QiZa) - F{Zs)^ with s e [m]), and thus: 

N m-2 . X 

+ lj =m!-(2m-l)(m-2)!-l- ^ft(^'^j(/i-l)! 
The independent parameters of (6.5) are therefore all together: 2m! — (2m — l)(m — 

m 

2)! + m + l-EftC)(^-l)!- 
1 

7. The simplest cases 

Let us see what in concrete terms happens calculating the formulas obtained 
in sections 5 and 6 for the first values of m = 2, 3, 4. 

A) m = 2. raj are of ra2ai which can therefore assume the values rio and roo- 
Furthermore C — D = ^, cqo = ''^oo and cio = rio and so, applying (5.7), we have: 



no! 



which, for (1.1), is equivalent to Dn- 

B) m = 3. Tq,^. are of r^gajai- As seen in section 6, the Fq — {r2io, ^'iio}) -^i = 
{''200}? F2 = {rioo}, -F3 = {''010} and the homologous Cq,^. are eliminated. Further- 
more q{Zi) = q{Z2) = qiZs) = rooo and cqoo = '^ooo and so, applying (6.5), we have 
that: 



(7.2) 



n n n ^ \ 

mn) = E/oE/E-oo(-i)^"-^°-^-°2/oi/( ; 

(rooo!)'(/ + 3rooo + 2)! ^ ^ nlrpop! /3rooo + / + 2^ 

cooo!(3rooo + 2)! ^ ^^f- ^ ' U \ f , 

fo+f+rooo=n 
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and we find (1.3) again. 

C) m = 4. Tocj are of ra^azoi^oci- The Fq = {r32io, r-22io, ''1210, niio, ''2110, '^siio}, 
Fi = {r32oo,''220o}, F2 = {r3ioo,riioo}, Fs = {rioio, ^2010}, F4 = {rono, ro2io} and 
the homologous Cq^. are ehminated. Furthermore q{Z[) = riooo + ''1200, Qi^i) ~ 
r2ioo + '^oioo, q{Z"') = rooio + ?'30io, ^(^2) = ^0200 + ^^1200, 9(^2 ) = ^2000 + 
^^2100, QiZs) = ?'3000+'^30io, and ^(Zi) = roooo+r3ooo+r20oo+ro20o, 5(^2) = roooo + 
'^3000+'^0100+'"l000, QiZs) = roooo+?'20oo+''iooo+''ooio, q{Z4) = roooo+''oioo+''ooio + 
ro20o. We besides have that: ciooo = riooo+''i20o-ci20o, coioo = 'roioo+''2ioo-C2ioo, 

COOIO = ''^0010 + ''3010 — C3010, C02OO = '''0200+ri200 — C120O! C200O = '''2000 +''2100— C2100, 
C3OOO ^''3000 + '■3010 — C3010, CoooO = '''0000 + ^"3000+''2Q0Q+''0200— C3000 — C2000— C0200, 

that:_i? = {riooo,ri2oo,r2ioo, ''0100, ''ooio, ''3010, ''0200, ''3000, ''2000, ''oooo} and that: 
C - D = {ci2oo, C2100, C3010} and, applying (6.3), we obtain: 



(7 3) K[4: n) — ^ ^ |^_]^^4n-/o-riooo-roioo-rooiO-r3000-r2000-ro200 

/o /s 'b. C-D 

Qfo2h+h+fs+uf „ „ ^V/2 + roooo + r3ooo + roioo + riooo)! 

V/O, Jl, J2, J3, J4, -K/ 

• if 3 + ''0000 + ''2000 + ''1000 + ''ooio)!(/4 + ''oooo + ^0100 + rooio + ''0200)! 

(/l + roOOO + ''3000 + ''2000 + ^0200)' 



coooo!c3ooo!c2ooo!co2oo! 



■('■0200 + '■1200)! (^2000 + ''2100)! 



/ , X ,Mooo + ''i20o\ /r2ioo + roioo\ /''oolo + '■30io\ 
('■3000 + '■3010)! 

V Ciooo,Ci200 / V C2lOO,CoiOO / V ^0010,^3010 / 



^ ^ ^ ^ ^_]^^/o+''1000+T'0100+'"0010+T'3000+T'2000+^0200 

fo+Y, fs+Y^R=ri C-D 

4 

«foo?''^' ^' , /''2000 + ^2100^ ./'■3000 + ^3010^ , 

6-'"2 1 Ci200'| _ _ C2l00!l _ _ )C3010! 



/o!n./Jn^!coooo! 
1 



V C2000,C2100 / V C3000,C3010 / 



_ / riQQQ + ri2oo\ / ''2100 + ''oioo\ / ''ooio + ''30io\ / ' 

\ Ciooo,Ci200 / V C2100,CoiOO / \ C0010,C3010 / \ 



''0200 + ''1200 
C02OO, C120O 

(/l + '' 0000 + '' 3000 + ''2000 + '"0200) K/2 + '' 0000 + '"3000 + '"OlOO + '' lOOo)' 
■ (/3 + ''0000 + ''2000 + ''1000 + ''001o)!(/4 + ''0000 + ''0100 + ''0010 + ''0200)! 



that is the result already obtained by Pranesachar and others in [1] . If instead we 
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apply (6.5) ,we obtain: 



(7.4) i^(4,n) 



fo+f+^R=n C-D 

6^°2f — =-(roooo + rsooo + ^2000 + ^-0200)! 




/0+''1000+''0100+T"0010+?'3000+'r2000+f0200 



"0! 




• ('^oooo + ^'aooo + ^0100 + ?^iooo)!(^oooo + ?'2000 + '^lOOO + '^ooio) 

■ (^ 0000 + ^0100 + '^OOIO + '^020o)K''lOOO + ^'l20o)K'"2100 + ^OlOo) 

• ('^ooio + »'30io)!(^02oo + ri2oo)!(r2ooo + r2ioo)!(r3ooo + rsoio)! 



/ + 4roooo + 2 (rsooo + r2ooo + roioo + riooo + ro2oo + rooio) + 3 



) 



Which is an improvement on the results known up to now, since it needs only 
15 independent variables (the ten of the three of C — D and the two /, /o) as 
compared with the 18 of the formula of Pranesachar, Athreya and Singhi. 



In conclusion we want to show how the Theorem 5.1 can have another inter- 
pretation which sheds light on its combinatory nature in a more profound way. 

The circumstance — which won't have escaped a careful reader — that the 
elements of R and of C are as many as those of &rn, and that is ml, is not casual. 
In fact if we interpret the indices am, ■ ■ ■ ,cxi and (3m, ■ ■ ■ , (^i as the inversion tables 
of one of the permutations of [m], putting bi = am+i-i (or bi ~ /3m+i-i), we will 
have two bijective maps between C and R and 6^, since < am+i-i < m — i. 

Furthermore C{c(j) will be the same as the number of cycles of cr e &m which 
corresponds in this way to Tq,^. . However it is not true — as could be thought — that 
the components Zi of aj correspond, in some way, to the cycles of the permutation 
a corresponding to Tq, . . 

To achieve this result we must introduce a new concept. Let us take a, a E &m 
written in its standard representation and put, Vi G [m] ki equal to n + l — t where t 
is the element furthest on the right among those to the left of i satisfying t > i {or if 
i = Sh, ki = m+l — St with st > and t maximum); moreover we set A;^ = if there 
are no elements greater than i on the left of i. We say that K{a) = (fci, . . . , km) is 
the "covering table" of a. It can be proved that the function K{a) is a bijection. 
Furthermore it is clear that < ki < m — i, Vie [m], and that, if fci = 0, i is a 
left-to-right maximum of the standard representation of a. 

Now, if we put ki = am+i-j, we have that. Vie [m], < ki < m — i and there- 
fore that (fci, . . . , km) can be interpreted as the covering table of a S{rctj) e &m- 
It can be easily proved that S{raj) is a bijection between R and &m and that, 
in this case too, C{aj) is the number of the cycles of S{raj)- Here however, if 
ShSh+i ■ ■ ■ Sh+p are the elements of a cycle of a written in its standard representation 



8. Another point of view 
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and if we take ks^,ks^+j^,. . . , kg^^^,, we have that {a^+i-s^, ccm+i-s^+i, ■ • ■ , am+i-sh+p} 
constitute a component Zm+i-Sh ^j- 

In the hght of this new bijective map, the results obtained previously can be 
expressed in a new combinatory language. In fact we can now consider the new 
variables To- and c^^ the indices of which consist of elements of &m (c, G &m) and 
again indicate their sets with R and C. Furthermore, writing 7 | o" G &rn to say that 
7 is a cycle of cr, we can put (5(7) = {r^ | 7 | cr} and 5(7) = ^ Q{l)'i corresponding 
meaning, going from to c^, will have Q(7) and ^'(7). In this way (5.6) can be 
reformulated like this: 



(8.1) K{m,n)= J], J], ^.^^.^.^ - ^(t) 



U^\a (-9(7)) i 



where W indicates the set of all the Va- in which a has an odd number of cycles. 

The simplifications of section 6, which conduct us to (6.5), can also be read 
more clearly now. In fact Fq consists of all the rcr in which a is made up of only 
one cycle of order m, while Fg, with s G [m], is formed of those in which a has 
a fixed point, made up of the element m + 1 — s, and a cycle of order m — 1 which 
permutes the other elements of [m] . 

Furthemore (8.1) reminds somehow the result attained by Gessel [3]. 

9. The latin squares 

When m — n, we find ourselves facing the Latin squares, much more famous 
than the Latin rectangles for their applications in various branches of mathematics. 

The number of n x n Latin squares is usually indicated by L{n). 

If we put m = n in (5.6) and in (8.1), and, abandoning the condition that the 
first row is in standard form, we multiply everything by n!, we obtain the following 
elegant result which allows us to count of the number of Latin squares of any order. 

9.1 Theorem. 

(9.1) 

n 

M^; 2^ n i^c Y\R\T\c\ ~ ^« ni^mc! 



In which the 2n\ parameters R and C, and their totals q{Zi) and 5(7) previously 
defined, appear. 
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